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ABSTRACT 

Recent developments of the Cascade-Exciton Model 
(CEM) of nuclear reactions are briefly described. These 
changes are motivated by new data on isotope produc- 
tion measured recently in "reverse kinematics" at GSI 
for interactions of 208 Pb and 238 TJ at 1 GeV/nucleon 
and 197 Au at 800 MeV/nucleon with liquid 1 H. This 
study leads us to CEM2k, which is a new version of 
the CEM code that is still under development. The 
increased accuracy and predictive power of the code 
CEM2k are shown by several examples. Further neces- 
sary work is outlined. 

INTRODUCTION 

The design of a hybrid reactor system driven with 
a high current accelerator requires information about 
residual nuclides that are produced by high energy pro- 
tons interacting in the target and in structural materi- 
als. It is both physically and economically impossible 
to measure all necessary data, which is why reliable 
models and codes are needed. A model with a good 
predictive power both for the spectra of emitted parti- 
cles and residual nuclide yields is the Cascade-Exciton 
Model (CEM) of nuclear reactions [1]. An improved 
version of the CEM is contained in the code CEM95 
[2], which is available free from the NEA/OECD, Paris. 
Following an increased interest in intermediate-energy 
nuclear data in relation to such projects as Accelera- 
tor Transmutation of Wastes (ATW) , Accelerator Pro- 
duction of Tritium (APT), and others, we developed 
a new version of the cascade-exciton model, CEM97 
[3,4]. CEM97 has a number of improved physics fea- 
tures, uses better elementary-particle cross sections for 
the cascade model, and due to some significant algorith- 
mic improvements is several times faster than CEM95. 
It has been incorporated into the recent transport code 
system MCNPX [5]. 

The recent GSI measurements performed using in- 
verse kinematics for interactions of 208 Pb [6,7] and 238 U 
[8] at 1 GeV/nucleon and 197 Au at 800 MeV/nucleon [9] 
with liquid -"-H provide a very rich set of cross sections 
for production of practically all possible isotopes from 
these reactions in a "pure" form, i.e., individual cross 



sections from a specific given bombarding isotope (or 
target isotope, when considering reactions in the usual 
kinematics, p + A). Such cross sections are much eas- 
ier to compare to models than the "camouflaged" data 
from 7-spectrometry measurements. These are often 
obtained only for a natural composition of isotopes in a 
target and are mainly for cumulative production, where 
measured cross sections contain contributions nor only 
from a direct production of a given isotope, but also 
from all its decay-chain precursors. Analysis of these 
new data has motivated us to further improve the CEM 
and to develop a preliminary version of a new code, 
(called CEM2k), which we describe below. 

RESULTS 

First, we discuss briefly the basis of the CEM and the 
main differences between the improved cascade-exciton 
model code CEM97 [3,4] and its precursor, CEM95 [2]. 
The CEM assumes that the reactions occur in three 
stages. The first stage is the IntraNuclear Cascade 
(INC) in which primary particles can be rescattered 
and produce secondary particles several times prior to 
absorption by, or escape from the nucleus. The ex- 
cited residual nucleus remaining after the emission of 
the cascade particles determines the particle-hole con- 
figuration that is the starting point for the second, pre- 
equilibrium stage of the reaction. The subsequent re- 
laxation of the nuclear excitation is treated in terms 
of the modified exciton model of pre-equilibrium decay 
followed by the equilibrium evaporative third stage of 
the reaction. Generally, all three components may con- 
tribute to experimentally measured particle spectra and 
distributions. 

An important ingredient of the CEM is the crite- 
rion for transition from the intranuclear cascade to 
the pre-equilibrium model. In conventional cascade- 
evaporation models (like ISABEL and Bertini's INC 
used in LAHET [10]) fast particles are traced down to 
some minimal energy, the cutoff energy T cut (or one 
compares the duration of the cascade stage of a reac- 
tion with a cutoff time, in the "timelike" INC models, 
like the Liege INC [11]) which is usually about 7-10 
MeV above the interior nuclear potential, below which 
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particles are considered to be absorbed by the nucleus. 
The CEM uses a different criterion to decide when a 
primary particle is considered to have left the cascade. 

An effective local optical absorptive potential 
W op t. m od.( r ) is defined from the local interaction cross 
section of the particle, including Pauli-blocking effects. 
This imaginary potential is compared to one defined by 
a phenomenological global optical model W op t.exp.( r )- 
We characterize the degree of similarity or difference of 
these imaginary potentials by the parameter 

V =| (W opt .mod. ~ W op t. exp.) /W op t. exp. \ ■ (1) 

When V increases above an empirically chosen value, 
the particle leaves the cascade, and is then considered to 
be an exciton. Both CEM95 and CEM97 use the fixed 
value V = 0.3. With this value, we find the cascade 
stage of the CEM is generally shorter than that in other 
cascade models. 

The transition from the preequilibrium stage of a re- 
action to the third (evaporation) stage occurs when the 
probability of nuclear transitions changing the number 
of excitons n with An = +2 becomes equal to the prob- 
ability of transitions in the opposite direction, with 
An = —2, i.e., when the exciton model predicts an 
equilibration has been established in the nucleus. 

The improved cascade-exciton model as realized in 
the code CEM97 differs from the CEM95 version by 
incorporating new and better approximations for the 
elementary cross sections used in the cascade, using 
more precise values for nuclear masses, Q-values, bind- 
ing and pairing energies, using corrected systematics 
for the level-density parameters, improving the approx- 
imation for the pion "binding energy", V^, adjusting 
the cross sections of pion absorption on quasideuteron 
pairs inside a nucleus, considering the effects on cascade 
particles of refractions and reflections from the nuclear 
potential, allowing for nuclear transparency of pions, 
and including the Pauli principle in the pre-equilibrium 
calculation. We also make a number of refinements in 
the calculation of the fission channel, described sepa- 
rately in [4]. In addition, we have improved many al- 
gorithms used in the Monte Carlo simulations in many 
subroutines, decreasing the computing time by up to a 
factor of 6 for heavy targets, which is very important 
when performing practical simulations with big trans- 
port codes like MCNPX. 

The authors of the GSI measurements performed a 
comparison of their data to several codes, including LA- 
HET [10], YIELDX [12], ISABELA (see details and ref- 
erences in [8]), and the Liege INC by Cugnon [11], and 
encountered serious problems: none of these codes were 
able to accurately describe their measurements; most of 
the calculated distributions of isotopes produced as a 
function of neutron number were shifted in the direc- 



tion of larger masses as compared to the data. While 
in some disagreement with the measurements, the Liege 
INC and ISABELA codes provide a better agreement 
with the data than LAHET and YIELDX do. Being 
aware of this situation with the GSI data, we decided 
to consider them ourselves, leading to the development 
of CEM2k. 

First, we calculated the 208 Pb GSI data [6] with the 
standard versions of CEM95 and CEM97. As an ex- 
ample, Fig. 1 shows our first results obtained with the 
CEM95 code. Let us note that so far all CEMxx codes 
simulate spallation only and do not calculate any pro- 
cesses following fission, such as production of fission 
fragments and the possible evaporation of particles from 
them. When, during a simulation of the compound 
stage of a reaction these codes encounter a fission event, 
they simply tabulate it (allowing us to calculate the fis- 
sion cross section and the fissility) and finish the cal- 
culation of this event without a subsequent calculation 
of fission fragments. Therefore, results from CEM95 
shown in Fig. 1 reflect the contribution to the total 
yields of the nuclides only from deep spallation pro- 
cesses (successive emission of particles from the target), 
but do not contain fission products. This is explicitly 
reflected in smaller calculated cross sections for light 
nuclides that are produced partially or mainly by fis- 
sion. We will not discuss them here. To be able to 
describe nuclide production in the fission region, these 
codes have to be extended by incorporating a model 
of high energy fission (e.g., in the transport code MC- 
NPX [5], where CEM97 is used, it is supplemented by 
the RAL fission model [13]). 

One can see that though CEM95 describes quite well 
production of several heavy isotopes near the target (we 
calculate p + 208 Pb; therefore, for us 208 Pb is a target 
not a projectile as in the GSI measurements), it does 
not reproduce correctly the cross sections for lighter iso- 
topes in the deep spallation region. The disagreement 
increases with increasing distance from the target, and 
all calculated curves are shifted to the heavy mass di- 
rection, just as was obtained by the authors of GSI 
measurements with all the codes they used. 



The results of the CEM97 code are very similar to 
those of CEM95 shown in Fig. 1. Our collaborator, 
Dr. V. F. Batyaev of ITEP, Moscow, performed an ex- 
tensive set of calculations of the same data with sev- 
eral more codes (HETC, LAHET with both ISABEL 
and Bertini options, CASCADE, CASCADE/INPE, 
INUCL, and YIELDX) and got very similar results [14]: 
all codes disagree with the data in the deep spallation 
region, the disagreement increases as one moves away 
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Isotopic distributions of the products in Pb-208+1GeV protons: GSI+CEM95 
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Product mass 

Figure 1: Comparison of isotope production cross sections from 208 Pb + p interaction at 1 GeV/nucleon measured 
at GSI [6,7] (filled circles) with the CEM95 calculations (open squares). 
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from the target, and all calculated curves are shifted in 
the heavy-mass direction. 

Similar disagreements with the GSI data obtained 
by different authors using different codes at different 
laboratories suggested to us that some general features 
of these reactions may be treated wrong by all codes 
and we started to modify our CEM97 code to attempt 
to understand the source of the disagreement. 

Fig. 2 shows some examples of steps we took in try- 
ing to determine which physical parameters might af- 
fect the discrepancies. First, according to Prokofiev's 
systematics [15], the value of the fission cross section 
for the reaction p (1 GcV) + 208 Pb should be w 100 
mb, i.e., only about 6% of the total reaction cross sec- 
tion of 1600 mb, as used by CEM97. In the past, if we 
were interested in describing only emitted particles and 
yields of the spallation products, we would not take 
into account fission processes in calculations of reac- 
tions with such small fission cross sections, to make 
the calculations faster. We would never be able pre- 
dict cross sections with an accuracy better than 6%. 
However, since many of the isotopes are produced with 
quite small cross sections, there is the possibility that 
the variation of fission barriers with neutron and pro- 
ton numbers might lead to an isotopic fractionation ef- 
fect. Calculations ignoring fission are shown both in 
Fig. 1 for all isotopes measured at GSI [6] (CEM95) 
and in Fig. 2a, for production of Tl, Ir, and Tm, cho- 
sen as examples (CEM95 and CEM97). As the general 
agreement of calculations with the data is not good, we 
decided to look at the influence of fission on the pro- 
duction of isotopes in the spallation region. In Fig. 2b, 
the solid lines show calculations with CEM97 including 
fission, while dashed lines show the results from CEM97 
(from Fig. 2a) not including fission. We see that though 
for most of the isotopes with high yields the two cal- 
culations almost coincide, the difference between the 
two calculations is somewhat larger for production of 
neutron-deficient isotopes of Tl (and to a lesser extent, 
Ir). We obtain similar results for other isotopes. This 
difference is large enough to lead us to the conclusion 
that if we hope to describe correctly production of all 
isotopes from a nuclear reaction, we have to take into 
account fission processes, even when the fission cross 
section itself is quite small. All the following results cal- 
culated by CEM97 are done including the fission mode, 
though this will not be mentioned explicitly in the text 
again. 



We know that a currently unsolved problem of 
CEM97, CEM95, and, to the best of our knowledge, of 
most of other available codes, is the correct prediction 



of the complex particle yields. Usually, all codes under- 
estimate production of 4 He, and to some extent, also of 
t and d at these incident energies. Such underestima- 
tion of complex-particle emission could affect the pro- 
duction of the final residual isotopes, so we decided to 
investigate the magnitude of this effect for our case. We 
are working now on better inverse cross sections used 
in calculating the probability for emission of complex 
particles and other improvements to the CEM hoping 
to address the problem of complex-particle yields, but 
for the moment we will estimate the effect of complex 
particle emission on residual nuclide production with a 
simple approximation. Our experience from analyzing 
other reactions involving similar energies and targets 
tells us that usually we underestimate the multiplici- 
ties of 4 He by a factor of 2, of t by a factor of 1.2, and 
of d by a factor of 1.1. In our model, these particles 
are produced only at the preequilibrium and evapora- 
tion stages of a reaction and their emission is governed 
by their widths. So, we multiply "by hand" the calcu- 
lated widths of 4 He by 2, of t by 1.2, and of d by 1.1 to 
estimate the effect. Our results are shown in Fig. 2c) 
as dashed curves compared with standard CEM97 re- 
sults from Fig. 2b shown as solid curves. Surprisingly, 
the two sets of calculations differ by little, so a simple 
enhancement of complex particle emission won't solve 
the problem of most residual nuclide yields. (Of course, 
for a specific isotope where a emission is an important 
mode of its production, the effect IS significant.) 

The standard CEM97 predicts the following multi- 
plicities for n, p, d, t, 3 He, and 4 He: 9.2409 ± 0.0025, 
1.9624 ± 0.0011, 0.5518 ± 0.0006, 0.1593 ± 0.0003, 
0.0682 ± 0.0002, and 0.1743 ± 0.0003, while if we en- 
large the widths of 4 He, t, and d as described above 
we get, respectively, 9.0886 ± 0.0025, 1.9198 ± 0.0011, 
0.5693 ± 0.0006, 0.1714 ± 0.0003, 0.0629 ± 0.0002, and 
0.02743 ± 0.0004. In other words, enlarging the widths 
of 4 He, t, and d does indeed predict more 4 He, t, and d 
as it should, but at the same time fewer n and p, so that 
the net effect on the residual nuclides is minimal. We 
conclude that although the problem of complex-particle 
yields needs to be solved independently of the GSI data, 
solving this problem will not affect significantly the pro- 
duction of most of the final residual nuclides, so we can 
try to find ways of describing the GSI data without first 
solving the problem of complex-particle yields. 

One of the things we previously noted was the fact 
that the calculated curves for almost all isotopes are 
shifted in the heavy mass direction and the shift in- 
creases as the atomic number decreases from that of 
the target. In other words this means that for a given 
final isotope (Z), all models predict emission of too few 
neutrons. Most of the neutrons are emitted at the evap- 
oration stage of a reaction. So, one way to increase the 
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Figure 2: Example of several steps from CEM97 to CEM2k in analyzing production of 81 T1, 77 Ir, and 69 Tm isotopes 
from 208 Pb interactions with 1 H at 1 GeV/nucleon measured at GSI [3] (see details in the text). 
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number of emitted neutrons would be to increase the 
evaporative part of a reaction. In our approach, this 
might be done in two different ways: the first would be 
to have a shorter preequilibrium stage, so that more ex- 
citation energy remains available for the following evap- 
oration; the second would be to have a longer cascade 
stage of a reaction, so that after the cascade, less ex- 
citon energy is available for the preequilibrium stage, 
so fewer energetic preequilibrium particles are emit- 
ted, leaving more excitation energy for the evaporative 
stage. 

Preequilibrium decay in our model involves both 
emission of particles from particle-like quasiparticle 
states, and the evolution of the number of excitons 
via equations for changes in exciton number of An = 
+2,0,-2 [1]. We make the transition to evaporation 
when the number of excitons exceeds a critical number 
which depends on the nuclear mass and excitation en- 
ergy. One easy way to shorten the preequilibrium stage 
of a reaction in CEM is to take into account only tran- 
sitions that increase the number of excitons, An = +2, 
i.e., the evolution of a nucleus only forward the com- 
pound nucleus, during the equilibration. We would thus 
not take into account possible transitions backward, de- 
creasing the number of excitons An = —2 and transi- 
tions without changing the number of excitons An = 0. 
In this case, the time of the equilibration will be shorter 
and fewer preequilibrium particles will be emitted, leav- 
ing more excitation energy for the evaporation. Such an 
approach is used by some exciton models, for instance, 
by the Multistage Pre-equilibrium Model used in LA- 
HET [10]. Calculations in this modification to CEM97 
are shown in Fig. 2d by dashed curves, compared with 
the standard CEM97 results shown by solid curves. We 
see that we get a slight improvement and a shift of the 
calculated curves in the right direction, but the effect 
is too small and doesn't solve the problem. 

In Fig. 2e, we show the results of the second method 
of increasing evaporation. To enlarge the cascade part 
of a reaction in CEM we have either to enlarge the pa- 
rameter V or to remove it completely and use a cutoff 
energy T cut , as do other INC models. Our calculations 
have shown that a reasonable increase of V doesn't solve 
the problem. The dashed curves in Fig. 2e show calcu- 
lations with CEM97 not using the parameter V, instead 
using a cutoff energy of T cut — 1 MeV, again compared 
with our standard CEM97 results shown by the solid 
curves. One can see a significant improvement of the 
agreement with the GSI data in this approach. Only a 
little more neutron emission is required to get just a per- 
fect agreement with the data. We choose to achieve this 
by applying in addition to this approach the condition 
of taking into account only transitions with An = +2, 
discussed above in reference to Fig. 2d. Using both 



these conditions leads to the results shown in Fig. 2f. 
We call this approach CEM2k and see that it describes 
the GSI data very well. 

The behavior of a compound nucleus in CEM is gov- 
erned by its mass number, A, charge, Z, excitation en- 
ergy, E* , and its angular momentum, L (we do not 
have in CEM any information about the deformation 
of a nucleus after the preequilibrium stage of a reac- 
tion). It is informative to look at these quantities for 
the compound nuclei remaining after the preequilibrium 
stage for the different approaches discussed above. Ta- 
ble 1 shows values of the mean A, Z, E* and L of the 
compound nuclei calculated using all the discussed ap- 
proaches. One can see that as a result of the transi- 
tion from CEM97 (plots a and b) to CEM2k (plot f), 
the mean mass and charge of the compound nuclei are 
almost the same, while their mean excitation energy 
increases by about 40 MeV. This is the main factor 
leading to a larger number of evaporated neutrons, and 
as a consequence, to a much better agreement with the 
GSI data. 

Table 1: Mean mass number A, charge Z, excitation 
energy E* (MeV), and angular momentum L (h) of the 
compound nuclei formed after the preequilibrium stage 
of reactions calculated by the different approaches as 
plotted in Fig. 2 



Method 


A 


Z 


E* 


L 


CEM97 


193.9 


78.2 


58.7 


24.6 


(plots a and b) 










Tq ► 2r Q 


193.4 


78.0 


58.3 


24.8 


(plot c) 










Only A„ = +2 


196.4 


78.6 


86.4 


24.1 


(plot d) 










No V 


191.0 


77.5 


65.0 


21.3 


(plot e) 










CEM2k 


193.9 


78.1 


97.5 


20.8 


(plot f) 











In short, CEM2k has a longer cascade stage, less 
preequilibrium emission, and a longer evaporation stage 
with a higher excitation energy, as compared to CEM97 
and CEM95. 

Fig. 3 shows a comparison of the CEM2k results with 
all the data (excluding those involving fission-fragment 
production) on isotope yields measured at GSI for the 
reaction 208 Pb + p. There is a very good agreement for 
almost all isotopes. 



A comparison of the CEM2k results with predictions 
by CEM95 [2], LAHET-ISABEL [10], LAHET-Bertini 
[10], CASCADE/INPE [17], CASCADE [18], YIELDX 
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Product isotopic distributions in * U0 Pb+1GeV: GSI+ITEP+CEM2k 
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Figure 3: Isotopic production cross sections of elements between Z=64 and 83 in the reaction of 208 Pb on hydrogen 
at 1 GeV/nucleon. Filled circles show GSI data [6], open stars are recent ITEP data measured by the 7-spectrometry 
method [16], while open squares show our CEM2k results. 
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[12], and INUCL [19] codes together with the GSI [6] 
and ITEP [16] data from our recent paper [16] is shown 
in Fig. 4. One can see that CEM2k agrees best with 
the data of the other codes tested here. 



Finding a good agreement of CEM2k with the isotope 
production yields, we wish to see how well it describes 
spectra of secondary particles in comparison with its 
precursor, CEM97. Fig. 4 shows examples of neutron 
spectra from interactions of protons with the same tar- 
get, 208 Pb at 0.8 Gev and 1.5 GeV (we do not know 
of measurements of spectra at 1 GeV, the energy of 
isotope-production yield data). We see that CEM2k de- 
scribes spectra of secondary neutrons certainly no worse 
than does CEM97, but possibly a little better at larger 
angles, in a good agreement with the data. So this pre- 
liminary version of an improved CEM code, CEM2k, 
describes both the GSI data from 208 Pb interactions 
with p at 1 GeV/nuclcon and the spectra of emitted 
neutrons from p+ 208 Pb at 0.8 and 1.5 GeV better than 
its precursor CEM97. 



Besides the 208 Pb data discussed above, reactions in- 
duced by 197 Au at 800 MeV/nucleon [9] and 238 U at 1 
GeV/nucleon [8] were measured recently at GSI. We 
use CEM2k as fixed from our analysis of the 208 Pb 
data [6,7] without further modifications to calculate the 
197 Au measurements [9] . Our CEM2k results are shown 
in Figs. 6 and 7 together with standard CEM97 pre- 
dictions and calculations by LAHET-Bertini [10] and 
YIELDX [12] codes from [9]. We see that just as in 
the case of the 208 Pb data, CEM2k agrees best with 
the 197 Au data in the spallation region compared to 
the other codes tested here. For the production of nu- 
clides lighter than Tb (Fig. 7), where fission-fragment 
formation begins to contribute, CEM2k underestimates 
the data (it doesn't include a model of fission-fragment 
production, as we discussed previously). 




T Lab (MeV) 



The 238 U GSI data are still preliminary (with rela- 
tive errors of 40%), published only as figures [8] and 
not available yet in tabulated form. These data are 
of a great interest to us as the reaction mechanisms 
in the case of p + 238 U interactions differ significantly 
from those of 208 Pb and 197 Au targets. Again, accord- 
ing to Prokofiev's systematics [15], for the reaction p 
(1 GeV) + 238 U the fission cross section is about 1390 
mb, which makes up about 80% of the total reaction 



Figure 5: Comparison of measured [20] double differ- 
ential cross sections of neutrons from 0.8 and 1.5 GeV 
protons on Pb with CEM2k and CEM97 calculations. 



8 



CEM2k LAHET-ISABEL CASCADE/INPE INUCL 

CEM95 LAHET-BERTINI CASCADE YIELDX 




150 160 170 180 190 200 210 

Product mass 



Figure 4: Isotopic mass distributions for independent production of Tm, Ir, and Tl isotopes from 1 GeV protons 
colliding with 208 Pb. Squares are ITEP measurements [16], while stars show GSI data obtained in reverse kinematics 
[6]. Results from different codes are marked as: CEM2k— our results, CEM95— Ref. [2], LAHET-ISABEL— Ref. [10], 
LAHET-Bertini— Ref. [10], CASCADE/INPE— Ref. [17], CASCADE— Ref. [18], INUCL— Ref. [19], and YIELDX— 
Ref. [12]. 
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Neutron number 

Figure 6: Isotopic distribution of the spallation products from the reaction 197 Au + p at 800 A McV from mercury 
to hafnium. Open circles are the GSI data [9], CEM2k (thick solid curves) and CEM97 (thick dashed curves) are our 
present calculations, LAHET-Bertini (thin solid curves) and YIELDX (thin dashed curves) are results of calculations 
from [9]. 
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Figure 7: The same as Fig. 6 but for products from lutetium to neodymium. 
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Neutron loss 

Figure 8: Comparison of the production of projectile 
isotopes in gold and uranium interactions with liquid 
hydrogen targets as measured at GSI [8] (symbols) and 
predicted by CEM2k (curves). 

cross section of 1788 mb as calculated by CEM2k. This 
means that fission is the main mode for this reaction 
and almost everything in a model prediction of isotope 
production should be determined by how well the model 
describes the neutron-to-fission width ratio, T n /Tf. So, 
the 238 U GSI data will be very useful to test and im- 
prove the treatment of the fission mode by any code, 
including CEM2k. We enlarged the figures from [8], to 
get (approximate!) numerical values of measured cross 
sections, and calculated this reaction with CEM2k. 

As an example, Fig. 8 shows a comparison of our 
CEM2k calculations with the GSI data [8] on produc- 
tion of uranium isotopes from p (1 GeV) interactions 
with 238 U and of gold isotopes from proton interactions 
with 197 Au. The agreement of CEM2k with 238 U data 
is not as good as for 197 Au, or as observed above for 
208 Pb at 1 GeV/nucleon and 197 Au at 800 MeV/n. This 
means we still have room for CEM2k improvement, es- 
pecially in the treatment of the fission mode. At the 
same time, Fig. 8 demonstrates how different are the re- 
action mechanisms in 197 Au and 238 U: while 21 isotopes 
of gold were measured from 197 Au, only 10 isotopes of 
uranium were measured from 238 U, and their yield de- 
creases exponentially as one moves into the neutron- 
deficient region, a direct result of the dominant role of 
the fission mode in uranium. 



FURTHER WORK 

Besides the changes to CEM97 mentioned above, we 



also made a number of other improvements and refine- 
ments, such as imposing momentum-energy conserva- 
tion for each simulated event (the Monte Carlo algo- 
rithm previously used in CEM provides momentum- 
energy conservation only statisticaly, on average, but 
not exactly for each simulated event); using real bind- 
ing energies for nucleons at the cascade stage of a reac- 
tion instead of the approximation of a constant separa- 
tion energy of 7 MeV used in previous versions of the 
CEM; using reduced masses of particles in the calcula- 
tion their emission widths instead of using the approx- 
imation of no-recoil used in previous versions; coales- 
cence of complex particles from fast cascade nucleons 
already outside the nucleus. These refinements are im- 
portant physically, and they slightly improve the agree- 
ment of calculations with the data, but not so dramat- 
ically as the changes discussed above. For this reason 
we do not show here figures with examples. In addition, 
CEM2k is still under development and will still change. 
We need to analyze a lot of more available data with 
it, especially at lower incident energies. We plan also 
to incorporate better inverse cross sections and to solve 
a problem of overestimation of fission cross sections at 
energies above 1 GeV, observed recently in our prelim- 
inary calculations of the p (2.6 GeV) + Hg reaction 
[21]. When we complete CEM2k to a reasonable level, 
we plan to incorporate it into LAHET and to replace 
the present CEM97 in MCNPX and the present CEM95 
in the MARS code system [22]. 
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